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A  Comparison  of  Free-Space  and  Fiber  Mixer  Performances 


1.  Introduction 

A  conventional  heterodyne  detection  system  employs  mixing  between  the  signal  field  and  the  local 
oscillator  (LO)  fields  in  free  space.  In  this  system  not  all  the  optical  power  collected  by  the  receiving 
aperture  contributes  to  the  power  in  the  intermediate  frequency  (IF)  signal.  Only  that  portion  of  the 
received  power  which  is  in  the  same  temporal  and  spatial  mode  as  the  LO  will  contribute  to  the  IF 
signal  [1], 

In  a  heterodyne  detection  system  employing  a  single  mode  fiber  mixer,  the  signal  and  LO  fields  are 
perfectly  matched,  so  all  the  signal  power  propagating  in  the  fiber  will  contribute  to  the  IF  signal. 
However,  not  all  the  power  focused  onto  the  end  of  a  single  mode  optical  fiber  will  couple  into  the 
fiber’s  single  propagating  mode. 

In  order  to  compare  a  heterodyne  system  employing  free  space  mixing  between  the  LO  and  the 
received  field  with  a  heterodyne  system  employing  mixing  in  a  fiber  coupler,  we  need  to  compare  the  IF 
carrier-to-noise  ratio  (CNR)  from  each  system  for  identical  system  design  parameters.  In  order  to  do 
this,  we  need  to  be  able  to  compute  what  fraction  of  the  received  optical  power  actually  couples  into 
the  fiber.  The  bulk  of  this  report,  and  Section  III  in  particular,  attacks  this  computation. 

The  numerical  results  here  indicate  that  for  most  reasonable  radar  optical  designs,  the  free  space 
mixer  will  do  as  well  or  better  than  the  fiber  mixer.  However,  the  fiber  mixer  can  perform  as  well  as 
the  free  space  mixer  when  the  receiver  optics  f-number  is  approximately  matched  to  the  f-number  of 
the  fiber.  This  will  often  result  in  a  reasonable  choice  for  the  receiver  optics  design. 
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U.  Free  Space  Mixer  Performance 


The  IF  signal  power  from  a  free  space  mixing  system  is  not  difficult  to  calculate.  All  the  equations 
necessary  to  make  the  computation  are  already  well  known  and  the  following  discussion  merely  collects 
those  results  into  a  form  suitable  to  solve  our  specific  problem. 

We  will  match  the  local  oscillator  (LO)  field  (referenced  to  the  transmitting  aperture)  to  the 
transmitted  field  and  define  the  fields  at  the  transmitting  aperture  as: 

UT(p,  =  2(0  £(P>.  U-LOfa  0  =  ^ To  e32™IFt  i(p)  (2  -  1) 


where:  P j.  is  the  transmitted  power, 

P lo  *s  the  LO  power, 

£(p)  is  the  transmitted  beam  shape  at  the  transmitter  aperture  plane,  normalized  to  unit  area 

f) 

(that  is,  J  J  K(P)T  =  1  where  Ap  is  the  area  of  the  transmitting  aperture), 

s(t)  is  the  transmitted  pulse  shape  normalized  to  unity  time  average  power  (indicated  by  the 
notation  <  >  ): 


T 

<  U(0r  >  =  y  j  dt  |  s( < )  |"  =  1  (where  T  is  the  pulse  width),  and 


i 'jp  is  the  IF  frequency  in  Hz. 

It  is  custom  .y  in  much  of  the  laser  radar  literature  to  write  the  detector  IF  current  as  a  target 
plane  integral  by  back  propagating  the  LO  field  to  the  target  [2].  Doing  this,  the  signal  component  of 
the  IF  current  out  of  the  detector  is  [3]: 


Ht)  =  £-0\\  A  d?  *  -  b  ILtfa t  -  ¥)  HU*'  0  (2  -  2) 

where:  q  is  the  electronic  charge, 

Tj  is  the  detector  quantum  efficiency, 
h  is  Planck’s  constant, 
vQ  is  the  optical  frequency, 

A  is  the  wavelength  of  the  light, 

L  is  the  range  to  the  target, 
c  is  the  velocity  of  light,  c  =  Xu0, 

At  is  the  target  area, 

T(p,  t)  is  the  target’s  random,  complex  reflectivity, 

Uf(p,  t  —  - £■)  is  the  far  field  pattern  of  the  transmitter  at  the  target  plane,  and 
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ty*c(7>,  t)  is  the  local  oscillator  (LO)  beam  pattern  back  propagated  to  the  target  plane. 

Using  the  Green's  function  for  free  space  propagation  and  making  the  far  Held  assumption 
Aj~  <  XL,  the  transmitter  and  back  propagated  LO  fields  at  the  target  plane  are  approximately1: 

Ut(p,t)  ^  7^  exp(  -  | p  p)  t) 

(2  -  3) 

^-0^^exp(-^|pP)ai^,0 
where  k  =  2ir/X  and  it  must,  of  course,  be  true  that: 

3m£'‘>  =  FFs<»V«?>)|z 

XL  (2-4) 

stio&o  =  -f  • 

JL 

Substituting  these  quantities  into  equation  (2  —  2),  we  can  write  the  IF  signal  component  as: 

=  (2-5) 

For  a  stationary,  purely  speckle  target  the  target  statistics  are  [3J: 

Ef&Pl )]  =  0,  E[m]  )T(p2))  =  0,  E[T(P  J  )T*(p2)  )  =  )6(p1  -p2),  (2-6) 

where  £[]  indicates  a  statistical  expectation  (average)  over  the  random  variable  or  process  and  ^Ts(p)  is 
related  to  the  reflectance  of  the  target2  at  coordinate  p-  If  the  target  is  range  unresolved, 
approximately  stationary  over  T,  and  uniform  over  the  area  illuminated  by  the  transmitted  beam  and 

j  9  i 

the  time  lags  (the  i  —  f  and  t  —  ^  effects)  are  negligible,  we  can  use  these  statistics  to  find  the  mean, 
time-average,  IF  signal  power: 


But  to  complete  this  calculation,  £()  needs  to  assume  a  specific  form. 

Let  the  transmitted  beam  shape,  £(/>),  be  the  usual  Gaussian  function  with  a  1/e  beam  power 

^In  this  report,  field  functions  in  a  script  font  indicate  the  Fourier  transform  of  the  respective  field  function  in  the 
Roman  italic  type.  That  is,  —  ^p{il7'(P'0)  |  _  where  the  notation  5p(L(j(P.  0)  |  _  means  take  the  Fourier 

transform  of  the  given  function  on  variable  p  and  evaluate  it  at  the  spatial  frequency  argument  w.  In  this  case. 
IL  j(p.  0  is  the  field  in  the  transmitter  aperture  plane.  See  Appendix  I  for  more  information  on  this  notation  and  how 
it  relates  to  other  notations. 

2For  a  uniform  target,  is  related  to  the  commonly  used  diffuse  target  reflectivity,  p,  by  ^  [9). 
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radius  of  w0,  normalized  to  unit  area  and  let  the  transmitting  aperture  A  y  be  large  enough  that  it 
does  not  appreciably  truncate  the  transmitted  beam  or  the  LO.  Then  the  the  beam  shape  and  its 
Fourier  transform  are: 


-'fr)  »  (2-8> 

W° '  XL 

If  the  target  is  resolved  by  the  transmitter  aperture  (»,e.,  the  illuminator  beam  extent  at  the  target  is 
less  than  the  target  area)  the  integral  over  the  target  is  easy  to  do  and  the  result  gives: 


|(P)  = 


ext 


irw‘ 


<E[\ 


PLOx2c*sPT 


(2-9) 


To  find  the  carrier-to-noise  ratio3  (CNR)  out  of  the  detector,  assume  the  post-detection  filter 
bandwidth,  B^,  is  greater  than  the  bandwidth  of  the  transmitted  waveform  s(t)  and  divide  by  the 

9 

expected  value  of  the  LO  shot  noise  component,  (q~r}/hv0)P ^qB ^  [1]: 


CNR 


free-space 


_ _ 2_ 


n  cr  nwo 
hv0BNPT**  [2 


(2-10) 


We  will  compare  this  with  the  CNR  from  a  fiber  mixing  system,  CNRjj^  ,  for  identical  system 
parameters  By,  Py,  Ts,  w0,  L,  and  Ay. 

It  is  interesting  to  examine  this  CNR  in  the  following  manner.  The  first  portion  of  the  CNR 

expression,  r///u/0B^ y,  is  just  the  inverse  of  the  IF  LO  shot  noise  power.  Therefore,  the  remainder  of 

2  o 

the  CNR  expression,  PyT^iru^/B  ,  is  the  IF  signal  power.  However,  the  expected  value  of  the  total 
power  collected  by  the  receiving  aperture  is: 

E[PC]  =  pT  ■  *- 7T  =  Vs~J-  (2-11) 

L  L 


Obviously,  only  a  portion  of  the  power  collected  by  the  receiver  aperture  actually  contributes  to  the  IF 
CNR.  Dividing  the  last  two  quantities,  assuming  a  monostatic  system  (so  A^  =  Ay),  and  using 
w0  =  dy/4,  (the  maximum  usable  value  of  w0  for  which  beam  truncation  is  negligible),  gives: 


1 

AR 


(2-12) 


In  other  words,  in  a  free-space  mixing  system  sensing  a  purely  speckle  target,  only  £  of  the  collected 
optical  power  actually  contributes  to  improving  the  IF  CNR.  The  remaining  collected  optical  power  is 


3 

There  is  a  difference  between  carrier- to-noise  ratio  (CNR)  and  signal- to-noise  ratio  (SNR).  The  CNR,  as  defined  here, 
is  sometimes  (and  incorrectly)  referred  to  as  the  SNR.  The  SNR  is  the  CNR  reduced  by  the  effect  of  variations  in  the 
received  power  due  to  target  variations  and  atmospheric  turbulence.  The  SNR  is  always  less  than  or  equal  to  the  CNR. 
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in  spatial  modes  which  do  not  mix  with  the  LO  spatial  mode.  It  is  possible  to  improve  on  the 
percentage  of  power  which  contributes  to  the  IF  CNR  by  using  a  Gaussian  beam  with  a  wQ  which  is  a 
larger  fraction  of  dj,  but  this  leads  to  beam  truncation,  increased  beam  divergence  in  the  far  field,  and 
the  above  equation  no  longer  r>  '  irately  predicts  the  coupling  efficiency  to  the  IF.  Rye  and  Frehlich 
have  shown,  however,  tha.  even  with  the  optimum  choice  of  u0  for  a  given  d j-,  the  maximum 
efficiency  when  using  truncated  Gaussian  beams  is  limited  to  about  44%  [10].  Similarly,  using 
dj/u’0  >  4  will  reduce  the  mixing  efficiency  to  even  less  than  1. 
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iH.  Fiber  Mixer  Performance 


A.  General  analysis 

Now  compute  the  IF  signal  power  for  the  fiber  mixing  system.  Since  the  local  oscillator  and  the 
received  optical  field  will  be  propagating  in  the  same  single  mode  fiber,  their  field  patterns  must  be 
identical.  A  single  mode  fiber  will  only  propagate  one  field  pattern,  called  the  HE^  mode.  Designate 
this  normalized  field  pattern  f/jj(/5).  The  signal  component  of  the  IF  current  out  of  the  detector  is: 


_  iv  ' 
h»0  \ 


PLOPf  *«-¥)« 


-  jllfUrpt 


IFt\  L 


P  p  s(t  _2I\ 
h^>\HLOFf  «*“-r  )e 


-  ib 


j2*vlFt 


(3-1) 


where  P j  is  the  received  optical  power  coupled  into  and  propagating  in  the  fiber,  Aj  is  the  detector 
area  (which  is  much  larger  than  the  mode  field  extent),  and,  once  again,  jZjj(p)  is  normalized  to  unit 
area.  Pj  is  a  random  variable  because  of  target  speckle.  Consequently,  the  mean,  time-average  IF 
signal  power  and  CNR  are: 


<£ll®‘)p)>  =(i^)2',toE[?/]  =>  CNRfibe,  =  (3'2) 

and  all  the  received  power  coupled  into  the  fiber  contributes  to  the  mean,  time-average  IF  signal 
power.  It  now  only  remains  to  compute  E[P  j]. 

To  compute  P j,  propagate  the  transmitted  field  U_F(p)  to  the  target  using  the  free  space  Green’s 
function, 


hL{p)  = 


(3-3) 


and  multiply  by  the  target’s  complex  reflectivity,  T(p).  Back  propagate  the  resulting  field  to  the 
receiving  aperture  using  the  same  Green’s  function  and  then  through  the  lens  to  its  focal  plane.  This 
will  give  us  the  field  at  the  lens’  focal  plane,  U_f(Pf),  in  terms  of  the  target  and  system  characteristics. 
The  end  of  the  fiber  will  be  at  the  lens’  focal  plane  where  it  will  collect  whatever  power  it  can  from  the 
field.  Finding  U_f(pj)  is  not  difficult.  Using  the  free  space  Green’s  function,  assuming  far  field 
operation  (Aj.  <  XL),  and  ignoring  insignificant  quadratic  phase  terms,  the  result  is  (see  Appendix  II 
for  a  complete  derivation  of  this  formula): 
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1 

(>A/)(AX)2 


(3-4) 


where:  /  is  the  focal  length  of  the  receiving  lens, 

— 7>(jZ’*)  is  the  Fourier  transform  of  the  beam  in  the  transmitter  aperture  (which  is 
proportional  to  the  far  field  beam  pattern  of  the  transmitter  at  the  target  plane),  and 

Wj^w)  is  the  Fourier  transform  of  the  receiving  aperture  function,  W j^(p).  (W pif)  is  1 
where  the  aperture  is  clear  and  0  where  it  is  opaque.) 

In  words,  this  formula  represents  an  image  of  the  target  in  the  focal  plane  of  the  lens,  with 
illumination  given  by  the  far  field  pattern  of  the  transmitted  beam,  and  with  the  image  blurred  by  the 
diffraction  introduced  by  receiving  aperture’s  Fourier  transform.  (This  is  easiest  to  see  by  letting  the 
receiving  aperture  be  infinitely  large  so  the  function  W ^()  becomes  a  Dirac  delta  function.) 

This  is  the  field  outside  the  fiber  face.  Inside  the  fiber  face,  the  field  is  that  of  the  HEjj  mode  of  a 
dielectric  waveguide.  The  field  distribution  of  this  mode  is  proportional  to  Jg(r)  in  the  fiber  core 
(where  r  is  the  radial  distance  from  the  center  of  the  fiber)  and  proportional  to  Kq(t-)  in  the  fiber 
cladding,  where  Jq()  is  the  zero  order  Bessel  function  and  Kq()  is  the  zero  order  modified  Bessel 
function.  Fortunately,  however,  Marcuse  (4], [5]  has  shown  that  a  Gaussian  function  very  accurately 
approximates  the  field  distribution.  Therefore,  the  normalized  single  mode  field  distribution  is 
approximately: 


LLU(P) 


(3-5) 


where  Marcuse  also  shows  the  optimum  choice  for  the  parameter  w  is  given  empirically  by  the 
equation: 


v  =  rcg0(V)  =  rc 


^0.65  + 


1.619  .  2-879  \ 
V3/2  V6  / 


(3-6) 


where:  rc  is  the  fiber  core  radius  and 

V  is  the  normalized  frequency  in  the  fiber,  V  =  £rc,'j2nAn~  (where  n  is  the  index  of  the  core 
material  and  An  is  the  difference  between  the  indices  of  the  core  and  cladding). 

Snyder  [6]  has  shown  that  it  is  possible  to  approximate  the  coupling  coefficient  between  a  field  in 
free-space  and  a  fiber  mode  with  an  overlap  integral.  Using  Snyder’s  general  formula,  the  power 
coupled  into  the  fiber  is: 
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(jA/)(AI)J 


l  \  ip  W)  31T(^)  } \<lf  !(*-)  VR(  -  f-jf) 


(jA/)(AI)2  J  J  At 


jj  &  mi  aT(i)  3t;,(  -  0)} 


P_ 
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where  the  second  equation  comes  from  writing  out  the  equation  for  U_j{)  and  rearranging  the  terms. 
The  third  equation  comes  from  writing  out  the  equation  for  W^()  and  rearranging  the  terms.  Defining 
the  truncated  or  apertured  version  of  the  fiber  field  mode  by  the  term: 


(3-8) 


XL 


and  using  it  in  the  previous  equation  gives: 


1 

(jXf)(XL)2 


(3-9) 


This  equation  has  a  very  satisfying  and  intuitively  appealing  form.  It  is  easy  to  see  how  the  form 
of  this  equation  is  very  similar  to  the  integral  portion  of  equation  (2  —  4)  for  'y(t)  for  the  free  space 
mixing  system.  The  first  term  here  is  the  target,  exactly  the  same  as  the  free  space  integral.  The 
second  term  is  the  transmitter  illumination  propagated  to  the  target  plane,  exactly  the  same  as  the  free 
space  integral.  In  the  free  space  integral,  the  remaining  term  is  the  local  oscillator,  referenced  to  the 
receiver  aperture,  back  propagated  to  the  target  plane.  Here  the  last  term  is  the  fiber  mode  field  back 
propagated  to  the  receiver  aperture,  truncated  by  the  receiver  aperture,  and  then  the  result  back 
propagated  to  the  target.  (The  Fourier  transforms  come  from  far  field  propagation  integrals.) 

Using  the  target  statistics  of  equation  (2  -  6)  and  taking  the  expectation  of  P  j  gives: 


(jXf)(XL)2 


The  last  term  in  the  integral  is: 

kht(&f  =  j  \*?\  WR&\)  exp(  f)‘ 

\\dp'2WRtf2)  atn(-0)exp(^|p-^f) 


(3-10) 
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where  the  approximation  is  good  for  the  far  field  and  the  change  of  variables  is: 


P0  -  1  2  "  ’  Ap  =  p'l  -  P2  =>  J  Jrfp'i  J  1^2  =  n^PoJ  JrfAp- 

Finally,  plugging  equation  (3  —  11)  into  equation  (3  -  10)  and  doing  some  algebra  gives  the  expression: 


El?/1=7W  I  ^  j 


A p 

XL 


WR&o  ~  \*p)  ^1 1(  ~ 


'  "W, 

-r-). 


(3-12) 


This  is  about  as  far  as  the  expression  can  go  without  giving  the  field  and  aperture  functions  some 
specific  forms. 


B.  General  circular  aperture 

Now  let  the  apertures  have  the  usual  circular  shape  with  diameter  dR: 

WRm  =  cir/fW'  m-dR'\ 

R  Vr)  \o  \p\>dRn 

let  the  transmitted  beam  be  the  usual  Gaussian  shape,  as  before: 


(3-13) 


£j(p)  — 


^  -Ti&=  nI 2™'°pt  exp(-© (3'14) 

and  the  fiber  mode  field  is  as  shown  in  equation  (3-5).  Plugging  all  this  into  equation  (3-12), 
simplifying,  and  scaling  the  space  variables  by  d ^  gives: 

||<^Poexp(-2(-xp)  |Pop]circ(Po  +  5Ap)circ(po-5Ap}  (3-15) 
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o 

Now  define  the  parameter  a  by  the  relation: 


-  2  CT/Rf 


(3-16) 


This  is  a  useful  parameter  because  it  appears  repeatedly  in  the  expression  for  E[P j]  and  it  is 
proportional  to  the  fiber  design  parameter  te/A  and  inversely  proportional  to  the  system  optics  f- 
number,  f/dft-  Now  the  expression  for  E[P j]  simplifies  to: 


/V  4T  qP T'd 


•  a2  J  J  dAp  exp(  -  ^  •  (^g)  +o2j|App). 
W^o  exp(  ~  4fl2  |  Po  f )  circ( Po  +  jSp)  circ(?o  “  gAp) 


(3-17) 


The  two  exponentials  are  just  standard  Gaussian  functions  centered  at  the  origins  of  the  pQ-  and 
Ap-planes.  The  circQ  functions,  however,  are  unit  diameter  circles  centered  at  ±  |a p  in  the  pG-plane 
and  the  p0-integral  is  just  over  the  area  of  overlap  of  the  two  circ()  functions.  If  |  Ap  |  >  1,  however, 
the  two  circ()  functions  will  not  overlap  at  all  and  the  entire  equation  is  identically  equal  to  zero. 
Figure  1  depicts  a  graph  of  circ^p0  +^Ap  j  and  circ^p0-iAp^  in  the  pc-plane  for  an  arbitrary  A p. 
Notice  the  area  of  overlap  will  always  be  centered  at  the  origin  of  the  pG-plane.  The  area  of  the 
overlap  is  dependent  only  on  (A p\  and  not  on  the  angle  of  Ap.  Furthermore,  since  the  Gaussian 
function  in  p0  is  circulo-symmetric,  the  inner  integral  over  the  p0-plane  will  also  be  independent  of  the 
angle  of  Ap. 

Figure  2  depicts  one  quadrant  of  the  p0- plane  for  the  specific  case  where  Ap  lies  on  the  y-axis, 
Ap  =  (o,  —  |  Ap  |),  so  the  center  of  the  circle  whose  arc  is  depicted  in  figure  2  is  at  coordinates 
(xc,yc)  =  ^0,  —  |Ap/2|j.  From  the  equation  for  a  circle  with  center  at  (xc,yc)  and  radius  r: 


(x  -  xc)2  +  (y  -  yc)2  =  r2, 


(3-18) 


it  is  easy  to  find  values  for  the  coordinates  on  the  x-  and  y-axes  of  the  arbitrary  point  on  the  edge  of 
the  overlap  area,  (x0,  y0),  shown  in  figure  2: 

xmax  =  5  A  1  i  Ap  P  J/maz  =  5  (*  -l^l)  =  2  ('i  *  —  4x°  _  I  O'  (3  —  19) 

Since  the  inner  integral  over  the  p0-plane  does  not  depend  on  the  direction  of  Ap,  it  is  possible  to 
simply  select  a  convenient  direction  for  Ap  and  perform  the  integration  for  that  direction.  Figure  2 
depicts  one  convenient  such  choice.  For  this  choice  of  Ap,  the  integral  over  the  entire  plane  is  just 
four  times  (since  there  are  four  quadrants)  the  integral  over  the  area  shown  in  the  figure. 
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Now  separate  the  integration  in  the  p0-plane  into  x-  and  y-integrals,  and  convert  the  outer  integral 
to  polar  coordinates  to  form  the  expression: 


(3  -  20) 

x  2 

where  erf()  is  the  standard  error  function  from  statistics,  erf(z)  =  \ Qdy  e-w  .  The  CNR  ratio  is 

then: 


There  are  two  parameters  in  this  expression,  dj^/w0  and  a.  Since  the  exp()  function  and  the  erf() 
function  are  both  very  smooth  functions  and  they  fall  off  rapidly  as  a  function  of  increasing  arguments, 
the  integrals  are  not  too  difficult  or  time  consuming  to  do  numerically.  Before  performing  this 
numerical  integration,  however,  let’s  examine  the  two  special  cases  of  very  large  and  very  small  a 
values. 

C.  Large  a  parameter  case 

If  we  let  the  a  parameter  be  large  enough  that  a  >  dR/u'o'  the  first  exponential  in  equation 
(3—17)  will  be  approximately  independent  of  d^/w0.  Furthermore,  if  a  is  this  large,  the  second 
exponential  will  be  sufficiently  concentrated  near  the  origin  of  the  p0-plane  that  the  p0-integral  will 
approximately  cover  the  entire  area  of  the  exponential  for  any  |  Ap|  <  1.  Then  we  can  approximate  the 
inner  integral  by  a  constant  (the  area  of  the  exponential)  when  |  Ap|  <  1  and  zero  elsewhere.  Therefore, 

E[Pj)\  ~4  1-T  fo-o2  J  \d5Tp  exp (  -  a2  [a^|~ )• 

J  1  large  a  irL“  J  J 

|  |dp0  exp(  -4o2  |pcf)  cir  <  Po  +  cir  c(p0-iAp) 
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- — fi  a2J  «xp( -o2  |i/>|")-  |  J <Jp0  exp( -4o2  |/>„p) 

=wi||^,xp(-02|^p) 


for  |  A p  |  <  1  and  zero  otherwise.  Now  perform  the  integral  over  the  unit  radius  circle,  |  A p  |  <  1  , 


_^sPTdR  n  1  —  e  ~  a  .  7sPTdR  x 
large  a  L2  2^  L 2  ? 


(3  -  23) 


where  the  approximation  assumes,  once  again,  a  large  a  so  that  e  a  <  1. 
Finally,  the  improvement  of  the  fiber  system  over  the  free  space  system  is: 


ttfiber(large  a) 
^Tree-space 


-f— f 

~\woa) 


(3-24) 


Unfortunately,  this  final  expression  is  greater  than  one  only  when  a  <  d ft/u>0  and  we  derived  this  result 
for  the  case  of  a  >  d^/w0\  Therefore,  in  the  large  a  parameter  limit,  the  fiber  mixer  will  always 
under-perform  the  free  space  mixer. 


D.  Small  a  parameter  case 

If  a<dft/w0,  the  first  exponential  in  equation  (3-17)  will  be  approximately  independent  of  a. 
Furthermore,  if  a  is  small  enough,  the  second  exponential  will  be  approximately  equal  to  one  and 
constant  over  the  region  of  the  p0-plane  enclosed  by  the  two  circQ  functions.  In  this  case,  the  inner 
integral  will  reduce  to  just  the  area  of  the  overlap  of  the  two  circ()  functions.  For  the  worst  case  of 
|  Ap |  =  0  (which  gives  the  greatest  extent  of  overlap  of  the  two  circ()  functions)  and  |p0|  =  |  this  will 
be  true  when: 

_2 

e  a  >  0.9  =»  a  <  0.3. 


The  overlap  area  is  just  twice  the  area  of  the  segment  of  a  circle  defined  by  an  arc  and  its  associated 
chord.  Figure  3  is  a  drawing  of  a  circle  »ment  with  the  arc,  chord,  and  other  parameters  labeled. 
For  our  problem,  r  =  1/2,  d  =|  Ap|/2,  and  h  =  r  —  d.  Using  standard  mensuration  formulas  for  circles 
(see  for  instance  [7]),  K  =  r^cos  —  ^(1  —  h/r)  -  (r  —  h)  N2rh  —  h?  and  the  area  of  overlap  in  terms  of  our 
variables  is: 

^(|A^|)  =  2A:  =  ^[cos-1(|A^t)-lA^|N|l-|  A^p  J  (3-25) 

Now  the  mean  value  of  the  power  propagating  in  the  fiber  is: 
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(3  -  26) 


~~2~^ '  *2  f  «*p(  -j(^)  l^f  )  -*(!  Ad) 


1  small  a  *L* 

Performing  a  change  of  variable  on  the  integral,  the  CNR  ratio  is: 


(3-27) 


Choosing  the  minimum  d^f wQ  value,  4,  and  performing  the  integration  numerically  gives  a  value  of 
~  0.0603  for  the  integral,  and 


CNRfiber(small  a) 

C^Rfree-space 


=  1.2285  a2. 


(3-28) 


Here,  once  again,  the  parameter  a  must  be  larger  than  about  1  in  order  for  the  fiber  mixer  to  out¬ 
perform  the  free  space  mixer.  But  we  derived  this  expression  under  the  small  a  assumption,  so  the 
fiber  mixer  will  under-perform  the  free  space  mixer  in  this  region  as  well.  However,  this  equation  and 
equation  (3  —  24)  together  indicate  the  possibility  that  there  is  an  intermediate  range  of  a  values  where 
the  fiber  mixer  could  perform  as  well  or  better  than  the  free  space  mixer.  The  only  way  to  determine 
what  happens  at  these  intermediate  f-numbers,  though,  is  to  perform  the  integration  of  equation 
(3  -  21)  numerically. 


E.  Numerical  analysis  results 

Figure  4  is  a  plot  of  the  CNR  ratio,  in  dB,  as  a  function  of  the  a  parameter  for  three  values  of 
dj^/ri’o,  4.0  (the  minimum  useful  value),  7.0,  and  10.0.  The  figure  includes  plots  of  the  asymptotic 
expressions  for  large  and  small  a  parameters  in  the  regions  where  they  are  appropriate.  Figure  4 
demonstrates  that  these  asymptotic  approximations  are  quite  good. 

Notice  the  fiber  system  is  always  worse  than  the  free  space  system  except  for  a  small  region  around 
an  a  parameter  of  about  two  where  the  improvement  is  on  the  order  of  a  few  dB  or  less.  At  least  part 
of  this  improvement  may  be  illusory,  however.  The  approximations  involved  in  forming  equations 
(3  -  5)  and  (3  -  7)  are  limited  in  their  accuracy.  In  fact,  Marcuse  [8],  in  plotting  his  figure  6,  found  a 
similar  excess  performance  peak  problem  when  using  exactly  these  same  two  approximations. 
However,  as  the  ratio  d^/ui0  gets  larger,  the  performance  improvement  becomes  more  pronounced. 
This  may  be  because  as  dj^/w0  gets  larger,  the  LO-received  field  matching,  as  given  by  equation 
(2  — 12)  gets  worse,  whereas  the  fiber  field-received  field  matching  may  remain  at  its  optimum  point 
with  the  appropriate  choice  of  parameter  a.  In  eflect,  the  choice  of  the  fiber  design  provides  an 
additional  parameter  to  vary  in  order  to  optimize  the  field  coupling.  Unfortunately,  this  is  more  a 
testimony  to  the  poor  coupling  for  large  values  of  d^/w0  than  a  praise  for  the  fiber  mixer.  In  practice 
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there  may  be  some  applications  where  a  large  d^/w0  ratio  would  be  useful.  An  example  might  be  a 
high  resolution  radar  (hence  the  large  d^)  which  used  flood  light  illumination  of  the  target  (hence  the 
small  wQ).  However,  for  cases  where  the  design  keeps  dR.I  w0  as  small  as  possible,  it  seems  the  free 
space  mixer  makes  the  best  possible  use  of  the  collected  backscattered  power  and  the  fiber  mixer  cannot 
significantly  improve  on  the  free  space  mixer’s  performance. 

Now  we  will  investigate  what  sort  of  fiber  and  optics  design  the  optimum  a  value  requires.  Recall 
the  a  parameter  is  defined  by: 


o 


a~ 


(3  -  29) 


and  the  ratio  w/\  is  defined  by: 


a,!*, 

A"  A  9°'V  1  ~  2ir^2nAn  2  x(NA) 


(3-30) 


where  N A  is  the  usual  numerical  aperture  for  a  fiber.  However,  for  a  single  mode  fiber,  the  range  of  V 
values  of  interest  for  is  just  1.5  <  V  <  2.4.  Within  this  range,  the  factor  Vg0(V)  is  numerically  almost 
constant.  It  varies  only  from  2.676  (for  V  =  1.5)  to  2.641  (for  V  =  2.4).  If  we  approximate  it  by 
Vg0{V)  ~  2.66,  we  find  w/X  ~  0.42/N  A  and  we  can  write: 


a  ~ _ 0^4 _ 

a-(NA)(f/m 


(3-31) 


where  //#  =  J  Mr  *s  ^e  standard  physical  v  jtics  f-number.  The  optical  f-number,  however,  is  just 
the  inverse  of  twice  the  numerical  aperture  [11].  Therefore  we  can  write: 


a  ~  1.9 


(^)optics 
(^fiber  ’ 


(3  -  32) 


From  figure  4  we  notice  the  optimum  chose  of  a  is  about  2  for  a  wide  range  of  d^/w0  values.  Our  last 
equation,  however,  tells  us  an  a  value  of  about  2  corresponds  to  the  point  where  the  numerical 
apertures  of  the  optics  and  the  fiber  are  approximately  matched.  This  is  exactly  what  we  would 
expect,  however,  from  simple  geometrical  optics  reasoning!  Furthermore,  since  reasonable  single  mode 
fibers  have  small  numerical  aperture,  we  will  require  receiver  optics  with  a  reasonably  large  f-number  in 
order  to  optimize  the  coupling.  This  is  an  fortuitous  result  since  the  target  the  optics  f-number,  the 
easier  it  is  to  fabricate  the  optics  with  low  aberrations. 

As  an  example  calculation,  return  to  equation  (3  —  31)  and  use  a  =  2  as  the  optimal  value.  Then 
we  get  the  best  performance  when 
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(3-33) 


//#  ~  2<JV'<W 

It  is  common  to  write  the  numerical  aperture  of  the  fiber  in  the  foilwing  form: 

^y‘)fiber  =  \lnl-n2-nl>f2A  (3-34) 

where  A  =  Reasonable  values  for  A  and  nj  (the  fiber  core  index  of  refraction)  for  a  single 

mode  fiber  are  A  =  0.003  and  =  1.5  which  gives  NA  ~  0.12  and  //#  2r  4.3.  Although  this  is  not  a 
large  f-number,  particularly  if  the  focal  length  must  be  long,  it  is  not  an  unreasonably  small  f-number 
either. 
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Figure  1 

CircQ  Function  Overlap 


eo-P^e. 
quadrant  I 


Figure  2 

CircQ  Function  Coordinates 


CNR 


Circle  Mensuration  [7] 


Figure  4 

CNR  ratio  -vs-  a  parameter 
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IV.  Summary  and  Conclusions 


This  report  compared  the  CNR  at  the  IF  stage  of  a  heterodyne  laser  radar  receiver  employing  the 
traditional  free  space  mixer  with  the  IF  CNR  from  a  heterodyne  laser  radar  receiver  employing  a  fiber 
mixer  stage.  The  analysis  shows  that  the  fiber  mixer  will  achieve  optimum  performance  when  the 
fiber’s  numerical  aperture  is  approximately  matched  to  the  receiving  optics  numerical  aperture.  This  is 
exactly  the  result  one  would  expect  from  simple  geometric  optics  reasoning.  When  the  fiber  and 
receiving  optics  are  optimized,  the  fiber  mixer  will  still  perform,  at  best,  only  marginally  better  than 
the  free  space  mixer,  however.  The  greatest  improvement  will  occur  in  cases  where  the  system 
requirements  force  an  inefficient  use  of  the  aperture  by  the  transmitting  beam.  In  this  case,  the  fiber 
mixer  can  compensate  somewhat  for  the  loss  of  mixing  efficiency  encountered  with  an  over  size 
aperture.  Apparently  the  free  space  mixer  makes  nearly  optimal  use  of  the  collected  optical  radiation 
even  though  the  mixing  efficiency  is  25%  or  less. 

Since  single  mode  fibers  always  have  quite  small  numerical  apertures,  the  optimal  receiving  optics 
in  a  fiber  mixer  may  have  a  reasonable  f-number.  A  typical  value  would  be  around  //#  ~  4.  This 
large  of  an  f-number  may  not  pose  too  much  of  a  problem  in  fabricating  low  aberration  optics.  It  does, 
however,  pose  an  additional  constraint  on  the  optical  system  design.  If  it  is  not  possible  to  use  the 
optimal  f-number  optics,  it  is  possible  to  use  the  relations  derived  in  this  report,  figure  4  in  particular, 
to  determine  how  much  of  a  reduction  in  mixing  efficiency  will  result.  This  will  give  a  better 
prediction  of  system  performance. 


-18- 


V.  References 


1.  Gagliardi,  Robert  M.  and  Sherman  Karp.  Optical  Communications,  John  Wiley  and  Sons,  New 
York,  1976. 

2.  Siegm&n,  A.  E.  “The  Antenna  Properties  of  Optical  Heterodyne  Receivers.”  Proceedings  of  the 
IEEE ,  54,  no.  10,  October  1966,  pp.  1350  -  1356. 

3.  Shapiro,  J.  H.,  B.  A.  Capron,  and  R.  C.  Harney.  “Image  and  target  detection  with  a  heterodyne- 
reception  optical  radar.”  Applied  Optics,  20,  no.  19,  1  October  1981,  pp.  3292  -  3313. 

4.  Marcuse,  D.  “Loss  Analysis  of  Single-Mode  Fiber  Splices."  The  Bell  System  Technical  Journal , 
56,  no.  5,  May  -  June  1977,  pp.  703  -  718. 

5.  Marcuse,  D.  “Gaussian  approximation  of  the  fundamental  modes  of  graded-index  fibers.”  Journal 
of  the  Optical  Society  of  America,  68,  no.  1,  January  1978,  pp.  103  -  109. 

6.  Snyder,  Allan  W.  “Excitation  and  Scattering  of  Modes  on  a  Dielectric  or  Optical  Fiber."  IEEE 
Transaction  on  Microwave  Theory  and  Techniques,  MTT-17,  no.  12,  December  1969,  pp.  1138  - 
1144. 

7.  Standard  Mathematical  Tables,  Sixteenth  Edition.  Samuel  M.  Selby,  ed.  The  Chemical  Rubber 
Company,  Cleveland,  OH,  1968. 

8.  Marcuse,  Dietrich.  “Excitation  of  the  Dominant  Mode  of  a  Round  Fiber  by  a  Gaussian  Beam.” 
Bell  System  Technical  Journal,  49,  October  1970,  pp  1695-1703. 

9.  Shapiro,  Jeffrey  H.  “Target  Reflectivity  Theory  for  Coherent  Laser  Radars."  Applied  Optics,  21, 
15  September  1982. 

10.  Rye,  Barry  J.  and  Rod  G.  Frehlich,  “The  truncated  Gaussian  lidar  antenna  problem  revisited,"  in 
Technical  Digest  on  Coherent  Laser  Radar:  Technology  and  Applications,  1991  (Optical  Society  of 
America,  Washington,  DC,  1991),  12,  pp  165-168. 


-19- 


Appendix  I 

Near  Field,  Far  Field,  and  Transform  Notations 


The  purpose  of  this  appendix  is  to  clarify  the  field  and  Fourier  transform  notations  used  in  this 
report  and  show  how  this  notation  relates  to  that  sometimes  used  in  other  literature,  like  that  in 
reference  [3]. 

If  we  define  the  near  field  distribution  at  the  transmitter  aperture  plane  by: 

UT(p,t)  =  fP^s(t)t(p)  (A1  —  1) 


where:  Pj.  is  the  transmitted  power, 

s(<)  is  the  transmitter  waveform,  normalized  so  <  |i(0p  >  =  1>  and 

£(p)  is  the  transmitter  beam  shape  normalized  so  that  J  J d~p  |£(/>)p  =  1  (where  Aj.  is 
the  area  of  the  transmitting  aperture), 
then  at  the  target  plane  the  far  field  distribution  is: 


(Al-2) 


where  C(p)  is  the  transmitter  beam  shape  at  the  target  plane  normalized  so  that  /  /  dp  |C(p)  f  =  1,  and 
is  given  by: 


C(p)=[[  df  i^hjjcp-p') 


(A1  —  3) 


-jkL  / 

where  k[Jj>)  is  the  free  space  Green’s  function,  k^ip)  =  ~J\l  exK  2X  1^1  r 

If  we  define  the  following  notation  to  represent  a  2-D  spatial  Fourier  transform  (as  in  the  footnote 
of  Section  II): 


(Al-4) 

XL 

we  can  represent  the  far  field  pattern  U f(p,  t)  in  terms  of  0  by: 

Ut(p,  t)  ~  ^  exp(  -  £  |  p  p)  3lT(^,  t)  ( A1  -  5) 

where  the  approximation  is  good  in  the  far  field  defined  by  Aj>  <C  XL. 

In  the  specific  case  of  a  Gaussian  transmitted  beam  and  a  transmitting  aperture  which  satisfies 
dj  >  4 wQ  so  the  transmitting  aperture  does  not  truncate  the  beam,  we  have: 
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£(?)  = 


<(p)  ~  jxf  exp( _  n  l *  ft  i2,r^  exp(-(!sr)2lpl2)“d 

—T^TL't)  -  -pT4^)*^2™©  «p(-©)  |Pp)>  (A1  —  6) 

where  once  again  the  approximation  is  good  in  the  far  field  defined  by  A  j  c  XL.  Here  w0  is  the  1  /e* 
beam  power  radius,  dj.  is  the  transmitting  aperture  diameter,  and  A  j  =  xdj/4  is  still  the  area  of  the 
transmitting  aperture. 
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Appendix  □ 

Derivation  of  Equation  3  —  4 


We  will  define  the  following  field  functions: 

U-tIPI'  0  **  Ike  transmitted  field  at  the  transmitter  aperture, 

J7f(7>2i  's  lke  transmitted  field  propagated  to  the  target  plane  (the  illumination  function  at 
the  target), 

Ur(P2 » 0  is  the  field  reflected  from  the  target, 

ILftiP 3*0  is  the  field  at  the  receiver  aperture  (tfr(/>2»0  back-propagated  to  the  receiver 
aperture),  and 

UApfJ)  is  the  field  at  the  focal  plane  of  the  lens. 

Furthermore,  the  free  space  Green’s  function  for  propagation  through  a  distance  L  is: 

ieJ*Iexp(?rl?P) 

Now  we  can  write  expressions  for  each  of  the  above  fields: 

0=||^  dP i  t  <  -  z)  -ifil  ~  ^2^’  (A^  -  1) 

Hr(P2'  0  =  U.0 2>  0  X  (P2*  0.  (A2  -  2) 

t)  =  1 1  ^  ^2  — r(^2’  *  ~  7)  £  -  l(?2  ~  ^3)’  ( A2  -  3) 

Uf(pf,t)=  |  dpzVR{p2,t)  ^7exp(^|?/|2).exp(-^/.p3)j  (A2  —  4) 

where  Ay,  Aj,  and  A^  are  the  transmitter  aperture,  target  area,  and  receiver  aperture,  respectively, 
and  /  is  the  focal  length  of  the  lens.  Notice  the  negative  range  value  in  the  Green’s  function  in  the 
equation  for  <)»  equation  A2  —  3.  This  is  because  this  expression  is  for  back-propagating  the 

reflected  field  to  the  receiver.  If  we  keep  the  same  coordinate  system  throughout,  the  field  is  now 
propagating  in  the  —  z-direction,  hence  the  negative  sign.  Plugging  one  expression  into  another  we  get: 

2/(?/.0=  {  \A  JH  [3^7  p/j2 -  j  ?/'ii3)]' 

J  J  ^  d?2  -  -  l^P 2  ~  3 )  X  (?2’  J  J  ^  ILt(p l» t  ~  %r)  ~ 

=  }  jjj  exp{57 1?/ f ) ' exp(  ~  J  ^/ ' ^3)  ‘ 
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1 1  a  d?2  *  -  L&2  ~  ^  (h> 1  -  T)  }  J  A  l  ~  ¥)  hL(P\  ~  ?2)-  (A2  -  5) 

where  W  ft(p3)  is  the  receiver  aperture  function  and  we  therefore  extend  the  limits  in  the  ^-integral  to 
the  entire  plane.  Writing  out  the  Green’s  functions  and  squaring  the  vectors  in  their  exponential  terms 
gives: 

=  w}(XL)t  •  ~<#/P>  1  \  ^3  "W  «p(  -jpfh} 

1  ~ (?2’ 1  -  i >  exp{  -  51  (l h  f  +  p3  f  -  2  *2  •  Pz))’ • 

rfpi£j(Pi,<-¥)exp(+|x(PlP  +  p2p“2^r^2))-  (A2  —  6) 

The  quadratic  terms  in  ^  cancel  while  the  quadratic  terms  in  pj  and  pg  are  insignificant  in  the  far 
field  where  Aj*  C  XL  and  A^  C  AL.  Furthermore,  the  quadratic  term  in  py  is  negligible  since  fiber 
core  diameters  are  on  the  order  of  5A  or  less  for  single  mode  fibers.  This  leaves  nothing  but  a  sequence 
of  Fourier  transform  integrals: 

Uf(Pf>  0  =  —-^2 11^3  Wr(p 3)  exp(  -  ^  Pf  •  p3)  J  J  ^dp2  2  (?2’  (-fo* 

cxp(t^2 '^3)}  #1  ^r(?l’*~¥)exp(-X?r^2)-  (A2-7) 

Under  the  assumption  that  we  design  the  system  such  that  the  diameter  of  the  transmitter  aperture  is 
large  enough  that  the  transmitted  beam  is  not  truncated  significantly,  we  can  extend  the  limits  of  the 
Pj-integral  to  the  entire  plane.  Then,  rearranging  terms  slightly,  we  have: 

=  I  A,02  2  CP2’‘  "MH  «p( -ihh} 

j  }  ^3  w  *^3)  exp{  -  7  Pf  Pz  )  •  «p(t  h  •  ?3  ) 

=  ^(^2  ll^i^P^-bll^Pi  VT(P\ •  *  -  ¥)  «p(  -  Wi  •  S)  • 

|  |  rfp3  exp(  -  j2wp3  •  — XX— ^), 

which  we  can  rewrite  in  our  Fourier  transform  notation  as: 
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2/0>/.0  =  (iV)(U)tf  f  Afh  1  '  -  M^I- '  -^)}k  •  *{»w}| 


XL 


L-,  - 

JPf  ~  P2 

XL 


(jXf)(XL)2i  J  At 

This  is  just  equation  3  -  4  of  the  text,  but  with  the  time  dependence  stated  explicitly. 


(A2  -  8) 
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